Phase space signatures of the Anderson transition 
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We use the inverse participation ratio based on the Husimi function to perform a phase space 
analysis of the Anderson model in one, two, and three dimensions. Important features of the 
quantum states remain observable in phase space in the large system size limit, while they would 
be lost in a real or momentum space description. From perturbative approaches in the limits of 
weak and strong disorder, we find that the appearance of a delocalization-localization transition is 
connected to the coupling, by a weak potential, of momentum eigenstates which are far apart in 
momentum space. This is consistent with recent results obtained for the Aubry-Andre model and 
provides a novel view on the metal-insulator transition. 



S3 
S3 

C/5 



I 

S3 
O 

o 



> 

CN 
CO 

o 

CO 

o 



i 

S3 
O 

o 



X 



PACS numbers: 05.60.Gg, 71.23.An, 05.45.Pq 



I. INTRODUCTION 

Phase space concepts are widely used in various areas 
of physics like quantum optics^ and quantum chaos^ 
while they are rarely employed in condensed matter 
physics. In this work, we use a phase space analysis 
to address the Anderson metal-insulator transition, and 
demonstrate that such a description is very useful and 
represents a powerful tool to describe and to elucidate 
how, as a function of a parameter, the nature of the 
eigenstates changes from delocalized to localized. 

While delocalized states call for a description in terms 
of momentum eigenstates, in particular in the ballistic 
regime, real space methods are expected to be appro- 
priate in the localized regime. Even though the real 
space wave function in itself contains already the full in- 
formation about a quantum state, a phase space repre- 
sentation may be much better suited to display the rele- 
vant information, e.g., in the vicinity of a delocalization- 
localization transition, where both, real space and mo- 
mentum space features, are expected to play an impor- 
tant role. 

The relevance of a phase space description has recently 
been illustrated by comparing the one-dimensional An- 
derson model and the Aubry-Andre model. In the first 
case, already the presence of very weak disorder leads 
to localized states in the thermodynamic limits In the 
quasiperiodic potential of the Aubry-Andre model, how- 
ever, a localization transition occurs at a critical poten- 
tial strength^ From a phase space analysis, it was con- 
cluded that this qualitative difference between the two 
one-dimensional models is due to the very different cou- 
plings of the momentum eigenstates, induced by the dis- 
order and the quasiperiodic potential, respectively^ 

In this work, we study the phase space behavior of the 
Anderson model in one, two and three dimensions and 
show that the above considerations are not restricted to 
one-dimensional models. In contrary, it turns out that 
the proposed relation between the coupling of momentum 
eigenstates due to a weak potential and the occurrence 



of a metal-insulator transition allows to explain why the 
Anderson transition cannot occur in one dimension. 

In view of the wealth of known results^ the Anderson 
model is particularly well suited for this kind of study. 
First studies of the one- and two-dimensional Anderson 
model based on the Wehrl entropy^ had already demon- 
strated that the diffusive regime present in two dimen- 
sions becomes apparent in phase spaced An extension to 
the three-dimensional Anderson model has become pos- 
sible by calculating inverse participation ratios in phase 
space instead of entropies. 11 Using the phase space analy- 
sis, we recover that in the thermodynamic limit all states 
are localized in one dimension while two dimensions rep- 
resent the marginal case. In three and higher dimensions, 
the phase space behavior provides clear signatures of the 
Anderson transition where states become localized only 
above a critical disorder strength. 12 This allows one to 
gain a detailed understanding of the phase space con- 
cepts and opens the road towards their application to 
more complicated systems. 

In Sec. ^ we start by introducing the characteriza- 
tion of quantum states by their inverse participation ra- 
tio (IPR) in phase space as well as the corresponding- 
quantities in real and momentum space. The Anderson 
model is introduced in Sec. IIIII and numerical results 
for its phase space behavior in the whole range from the 
ballistic to the localized regime are presented in Sec. IIVI 
The observed features are discussed in the light of known 
properties of the eigenstates. Since the behavior in the 
limiting cases of weak and strong disorder turns out to 
depend on the dimensionality and to be indicative of the 
existence of a metal-insulator transition, we devote the 
main part of this paper to a detailed investigation of these 
limits. Perturbative expansions for the inverse participa- 
tion ratios in the different spaces are presented for strong 
disorder in Sec. IV Al and for the limit of weak disorder in 
Sec. IV Bl Here, a crucial dependence of the inverse par- 
ticipation ratio on dimension is identified, and related to 
the structure of the coupling of momentum eigenstates 
by weak disorder. This important property is only ap- 
parent in phase space while such signatures cannot be 
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extracted from the inverse participation ratios neither in 
real nor momentum space. Our interpretation and the 
relation to the known properties of the Anderson model, 
in particular in the marginal case of two dimensions, is 
confirmed by an analysis of the dependence of the inverse 
participation ratio on system size in Sec. IVII Finally, we 
present our conclusions in Sec. IVIll 



II. CHARACTERIZATION OF STATES 

Among the infinite variety of possible phase space 
representations of a quantum stated the Husimii^ or 
Q function^ is best suited for our purpose because it 
guarantees a positive definite density. This property will 
allow us to define an inverse participation ratio in Eq. © 
below. The positivity is a direct consequence of the def- 
inition of the Husimi function 



p(x ,ko) = \(x ,k \ip)\ 7 



(1) 



where the state \ip) is projected onto a minimal uncer- 
tainty state |xo, ko) centered around position xo and mo- 
mentum ko in phase space. The minimal uncertainty 
state assumes a Gaussian form both in position and mo- 
mentum representation. Its real space wave function 
reads 



<*l x o,ko)=(^) 
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In the definition of the Husimi function, the width 
a appearing in J5J) determines the relative importance of 
structures in real and momentum space. We adopt this 
definition for lattice models with periodic boundary con- 
ditions provided that er <C L. Here, L is the number of 
lattice sites in one spatial direction and the lattice con- 
stant sets the unit length. Throughout this paper, we 
choose a = ^/L/Att which yields an equal width of the 
Gaussian relative to the system size L and the momen- 
tum interval running from k = — n to n. 

Since we are ultimately interested in the thermody- 
namic limit, L — > oo, let us first discuss the dependence 
on system size of the phase space resolution provided by 
the Husimi function. Since the d spatial components are 
independent of each other, it is sufficient to consider the 
one-dimensional case. For our choice of a, the Gaussian 
smearing arising from the projection onto a minimal un- 
certainty state affects areas in phase space which contain 
of the order of \f~L x \f~L grid points. Structures appear- 
ing on smaller scales cannot be resolved. However, even 
though the absolute resolution degrades, relative to the 
size of the system the resolution becomes increasingly 
better as the system size is increased. This holds for any 
(7 which scales with system size like L a where < a < 1. 
In contrast, the limiting cases a — and a = 1 behave 
quite differently. For a — 0, we have optimal resolution in 
real space but cannot resolve phenomena in momentum 
space, even in the thermodynamic limit. The opposite is 



true for a — 1 where one would obtain a pure momentum 
space description. Our choice of a = 1/2 leads to an ideal 
balance between these two extreme cases, and allows to 
track features which rely on both, real and momentum 
space, to the thermodynamic limit. 

The Husimi function contains a tremendous amount of 
information about a quantum state. It turns out, how- 
ever, that relevant information can already be extracted 
by considering the inverse participation ratio (IPR) in 
phase spaceii 



x.k 



(3) 



where the sum runs over all phase space points (x, k) . 
The normalization in is chosen in such a way that 
P = 1 corresponds to an optimal localization around 
one lattice point. In phase space, this is achieved by a 
minimal uncertainty state. A distribution of the Husimi 
density over a larger volume in phase space corresponds 
to lower values of P. 

Although the IPR in phase space J3J is defined in 
terms of the Husimi function p(x, k), it may be calcu- 
lated directly from the wave functionjii*i£*il by means of 
a straightforward genera lization of the one-dimensional 
expression given in Refs. Illlld Such an approach pro- 
vides significant numerical advantages and is crucial for 
the treatment of higher-dimensional systems. 

The IPR in phase space P should be compared with the 
IPR in real space which has frequently been employed to 
describe quantum states in disordered systemsii&iiLSSiSi 
Here, the state \ip) is projected onto a Wannier state |x) 
localized on a single site of the lattice. This allows to 
define the IPR in real space as 



Pr 



Ei 



(4) 



which corresponds to the limit a — > of the IPR in phase 
space. It is also convenient to introduce the IPR in mo- 
mentum space as 



p k =^2\(m 



(5) 



where the basis of momentum eigenstates |k) is given by 
(x|k) = exp(ik-x)/L d / 2 . 

As will be seen below, even the combined information 
from the IPRs in real and momentum space is not equiv- 
alent to the information provided by the IPR in phase 
space. However, it was shown in Ref. |22| that by an 
appropriate Gaussian smearing of the real and momen- 
tum space densities one can define marginal distributions 
which allow to reproduce the behavior of the IPR in phase 
space. Unfortunately, this approach does not result in a 
reduction of the numerical effort as compared to the cal- 
culation of the IPR in phase space. 
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III. THE ANDERSON MODEL 

In the following, we shall present a detailed compari- 
son of the IPRs in real, momentum and phase space by 
considering the Anderson model for a quantum particle 
in a disordered potential. Its Hamiltonian 

H = -t J2 (|x')(x| + |x)(x'|) + W/5>„|x)(x| (6) 

<x,x'> n 

is defined on a rf-dimensional square lattice with L sites 
in each direction. The energy scale is set by the hopping 
matrix elements t = 1 between nearest neighbor sites 
< x, x' >. In order to avoid boundary effects we choose 
periodic boundary conditions in each direction so that 
every site has 2d nearest neighbors. The on-site energies 
v n forming the disordered potential are drawn indepen- 
dently from a box distribution on the interval [—1/2; 1/2] 
and W denotes the disorder strength. 

The structure of the quantum eigenstates of the An- 
derson model depends on the disorder strength. For van- 
ishing disorder, the eigenstates are plane waves and thus 
are localized in momentum space. In the opposite limit 
of strong disorder, localization in real space takes place. 
In order to describe the behavior of the states in the 
whole parameter region, and in particular the transition 
between the limiting regimes, it is very useful to work 
with phase space quantities which adequately take into 
account real space as well as momentum space properties 
at the same time. 



IV. INVERSE PARTICIPATION RATIOS FOR 
THE ANDERSON MODEL 

In order to appreciate the advantage of the phase space 
approach, we start by comparing the IPR in real space, 
P x , phase space, P, and momentum space, P&, for the 
two-dimensional Anderson model. In Fig. ^ numeri- 
cal results are shown for a lattice of size 64 x 64. For 
each given disorder strength W, we have diagonalized 
the Hamiltonian © for 50 different disorder realizations 
{v n }, and used L 2 /2 states around the band center to 
calculate distributions of logarithms of the IPRs. 

In Fig. \Ql we observe a monotonic increase of the real 
space IPR with increasing disorder strength W. This 
corresponds to the tendency towards localization of the 
eigenf unctions. According to Fig. the IPR in mo- 
mentum space simultaneously decreases, thereby indicat- 
ing derealization in momentum space. This behavior of 
P x and Pk is an immediate consequence of the system's 
change from the ballistic regime for weak disorder, e.g., 
localization in momentum space, to localized states in 
real space for strong disorder. 

Since the IPRs in real and momentum space evolve in 
opposite directions as a function of the disorder strength, 
the behavior of the phase space IPR, which describes the 
spread of the wave function in real and momentum space 
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FIG. 1: Gray scale plot for the distributions of the logarithms 
of the inverse participation ratios in (a) real space, (b) phase 
space, and (c) momentum space as a function of the disor- 
der strength W. The data represent L 2 /2 states around the 
band center, for 50 different disorder realizations of the two- 
dimensional Anderson model of size L = 64. 



on an equal footing, can be expected to provide more 
subtle informations. Indeed, the behavior of the phase 
space IPR depends on the details of the model as can 
be seen by a comparison of the one-dimensional Ander- 
son model and the Aubry- Andre model^ and within the 
Anderson model itself, where the dimensionality plays a 
crucial roleiii 

For the two-dimensional case, the IPR in phase space 
depicted in Fig. ^> displays a much richer structure than 
the IPRs in real and momentum space. In particu- 
lar, the dependence on the disorder strength W is non- 
monotonic, and one finds a minimum at an intermedi- 
ate value of W which can be associated with diffusive 
behavior This non-trivial behavior motivates the fol- 
lowing in-depth study of the Anderson model by means 
of the IPR in phase space. 

Figs. |2Jl-c depict the mean IPR in phase space for one, 
two, and three dimensions, respectively, for various sys- 
tem sizes L. The color changing from red to black cor- 
responds to an increasing system size which is also indi- 
cated by the arrows. The data have been scaled with the 
length dependence L~ d / 2 of the limiting cases at W = 
and W — > oo, cf. Sec. El Before giving a detailed discus- 
sion of the dependence on L in Sec. lVfl we concentrate on 
the overall behavior as a function of the disorder strength. 

One of the most striking aspects of the results pre- 
sented in Figs. [2^.-0 is that the behavior of the phase space 
IPR at weak disorder depends on the spatial dimension 
in a crucial way. While in d = 1 the IPR increases with 
increasing W, it decreases in d > 2. Together with the 
fact that, independently of the dimension d, at strong dis- 
order the limiting value for W — ► oo is approached from 
above, this has important consequences for the global be- 
havior of the phase space IPR. In d = 1, the two limits are 
joined by a peak indicating localization in phase space. 
In contrast, in two and three dimensions, P decreases in 
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FIG. 2: Mean IPR in phase space as a function of the disorder 
strength for (a) one, (b) two and (c) three dimensions. The 
system size L increases from the red to the black curves: (a) 
L = 128, 192, 256, 384, 512, 768, 1024, 1536, and 2048; (b) 
L = 16, 24, 32, 48, 64, and 96; (c) L = 14, 16, 18, 20, 22, and 
24. The arrows indicate how the position W of the extrema 
shifts with increasing L. 



FIG. 3: The standard deviation a(P) of the IPR in phase 
space as a function of the disorder strength for (a) one, (b) 
two and (c) three dimensions and the same parameters as in 
Fig. El 



PERTURBATION THEORY 



the regime of small disorder, and assumes a minimum 
indicating a large spreading in phase space followed by 
a more or less steep rise towards a maximum as can be 
seen in Figs. Ep and c. 

The minimum of the phase space IPR in two and higher 
dimensions can be associated with the existence of a dif- 
fusive regime where the system size is much larger than 
the mean free path but smaller than the localization 
length. The resulting mixing of the plane waves by the 
disorder potential considerably alters the structure of the 
states and leads to a spreading both in real and momen- 
tum space and thus to a small value of the phase space 
IPR. This is reminiscent of the emergence of quantum 
chaos and can be confirmed by determining the energy 
level statistics around the minimum of P. One indeed 
finds the Wigner-Dyson statistics^ which characterizes 
the diffusive (chaotic) regime. 

In addition to the mean value, the distribution of the 
phase space IPR at given disorder strength (cf. Fig. 
for the case d = 2) can be characterized by the standard 
deviation a(P) depicted in Fig. for one, two, and three 
dimensions. Here, we have employed the same scaling 
with system size as in Fig. |5J The overall structure re- 
sembles the one found for the mean values. The strong 
suppression of the standard deviation occurring in the 
diffusive regime, being particularly pronounced in d = 3, 
indicates that the phase space structure is quite indepen- 
dent of the individual states. This confirms once more the 
universal chaotic character of the diffusive states which is 
expected due to the strong mixing present in this regime. 



The numerical results for the phase space IPR pre- 
sented in Figs. Eland indicate that the changes in the 
global behavior as a function of the disorder strength can 
be understood in terms of the limiting behavior for strong 
and, in particular, for weak disorder. Therefore, we pro- 
ceed next to a detailed perturbative investigation of the 
IPRs in these two limits. We start with the simpler case 
given by the limit of strong disorder. 



A. IPR at strong disorder 

For W — > co, all eigenstates are localized on single sites 
in real space. A finite ratio t/W then leads to a coupling 
to the nearest neighbor sites due to the kinetic energy 
in (JBJ. Such a perturbation can be treated analytically 
as long as it is sufficient to take into account only the 
coupling to the nearest neighbor state which is closest 
in energy to the initial site. For the resulting two state 
system, the IPRs may be calculated explicitly. The other 
nearest neighbor sites enter in the calculation only when 
the disorder average is performed. 

In a first step, we thus focus on two nearest neigh- 
bor Wannier states on a lattice of size L d . The absolute 
value of the difference between the corresponding on-site 
energies will be denoted by A. Then the effective Hamil- 
tonian for the two level system in the Wannier basis reads 



Htls 



— A/2 -t 
-t A/2 



(7) 



It is straightforward to determine the two eigenstates 
and the corresponding IPRs, which are identical for both 
states. Introducing the eigenenergies A = ±[(A/2) 2 + 
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i 2 ] 1 / 2 , the IPRs are given by 
t 2 

P*(A) = 1 



2A 2 



P fc (A) = L- d 
P(A) = iT d/2 



1 



2A 2 



(8) 



1 + -^(2exp(-l/4cr 2 ) - 1) 
2A 2 



In particular, one finds P x (0) = 1/2 because for degen- 
erate on-site potentials the two states are both equally 
distributed over the two sites. Furthermore, and consis- 
tent with the results of the previous section, the IPRs in 
real and momentum space behave oppositely as t/W is 
increased. For large system size, the IPR in phase space 
increases with t/W , just as the IPR in momentum space. 

In order to compare with our numerical results, we 
need to perform a disorder average. Since the on- 
site energies are equally distributed inside the interval 
[— W/2; W/2], the probability density pi that two neigh- 
boring on-site energies differ by A reads 



Pi(A) = w? (W-A). 



(9) 



The index 1 indicates that only one nearest neighbor site 
is taken into account. 

Furthermore, we need to ensure that the energy dif- 
ference A is the smallest among the energy differences 
with all nearest neighbors. Therefore, for the remaining 
2c? — 1 nearest neighbors, the difference in on-site energy 
with respect to the central site should be larger than A. 
The probability density for such a 2d nearest neighbor 
configuration is given by 



p 2d (A) = AT 1 Pl (A) 




2d-l 



dxpi(x) 



(10) 



where N is a normalization constant. 

Within the assumption that we can restrict ourselves 
to an effective two level system, we therewith obtain the 
IPR in phase space 



P = 



w 



dA P2d (A)P(A) 



(11) 



together with corresponding expressions for the IPRs in 
real and momentum space. Making use of © and IjlOfl. 
to leading order in t/W, one obtains 



Pk = L 



1 + 2nd± 



(12) 



P = L 



l + 2nd— [2exp(-l/4a 2 ) - 1] 
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FIG. 4: Comparison of the perturbative result 2nd(t/W) 
(solid line) with the numerically computed values for devia- 
tions of the IPRs from their value at W = oo in real space (c x , 
dotted line), momentum space (ct, dashed line), and phase 
space (c, dashed-dotted line), for the parameters of Fig. Q 



with corrections of order (t/W) 2 \n(t/W). Configura- 
tions where more than one nearest neighbor site is ener- 
getically degenerate with the central site do not modify 
the results Ijl2|l because the probability to find such a 
configuration vanishes. 

In order to compare the numerical data presented in 
Fig. ^ with the perturbative result, we introduce the 
quantities c x = 1 — P x , Ck = L d Pk — 1, and c — 
(L d l 2 P - l)/[2exp(-l/4cr 2 ) - 1]. Within the perturba- 
tive results of Eq. we have c x = Ck = c = 2ird(t/W). 
The numerical results for the two-dimensional Anderson 
model are shown in Fig.0] The agreement with the lead- 
ing perturbative results is remarkably good for disorder 
strengths down to rather small values of W. This is 
particularly true for the IPR in momentum space rep- 
resented by the dashed line. The fact that IPRs are by 
definition positive quantities implies that the correction 
c x of the IPR in real space depicted by the dotted line 
is limited from above by 1. Therefore, the leading cor- 
rection to P x given by Ijl2|l must fail when c x reaches 
this limiting value. Finally, the dashed-dotted line cor- 
responding to the phase space term c is well described 
by the leading perturbative correction according to l|12(l 
down to W ss 30 for this system of 64 x 64 sites. 

As for the case of two sites, Eq. (|SJ|, the IPRs in real 
and momentum space move in the opposite direction as a 
function of t/W. Moreover, the IPR in phase space still 
behaves similarly to the momentum space IPR, for the 
averaged quantities given by (|12|l . The key to an under- 
standing of this behavior of the IPR in phase space lies in 
the limited resolution provided by the Husimi function. 
Since its spatial resolution is of order L 1 / 2 , changes which 
occur only on two lattice sites will not affect the Husimi 
functions, in particular in the case of large system sizes. 
Only the small deviation of the factor 2 exp(— 1/4ct 2 ) — 1 
from 1 can be traced back to real space behavior as an 
incomplete overlap of the Gaussians centered at the two 
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sites in question. Small scale changes in real space, how- 
ever, lead to large scale changes in momentum space. 
In the regime discussed above, one observes beatings in 
the momentum space density as a consequence of the re- 
quired orthogonality of the two eigenstates of J7J. This 
effect can be resolved by the Husimi function, so that mo- 
mentum space effects dominate the phase space behavior 
at strong disorder. 

Finally, the difference of L d t 2 in the prefactor of the 
IPRs in phase space and momentum space stems from 
the Gaussian smearing in phase space which contributes, 
in our case of spatially well-localized states, a factor a oc 
L 1 / 2 for each spatial dimension. 



B. IPR at weak disorder 

As shown in Fig. ^ the IPRs in real and momentum 
space exchange their qualitative role as compared to the 
strong disorder limit. This is not surprising, because bal- 
listic motion of a quantum particle implies the existence 
of plane waves with well localized momentum and de- 
localization in real space. Exchanging real and momen- 
tum space, this corresponds to the real space scenario for 
strong disorder. 

The situation, however, is more complicated in phase 
space, and the behavior in the limit of weak disorder, 
W — ► 0, is by far more complex. Only in one dimen- 
sion, the IPR in phase space can indeed be understood 
in terms of the real space IPR at weak disorder and the 
momentum space IPR at strong disorder >Ai In particular, 
the IPR in phase space increases with increasing disorder 
strength in the regime of weak disorder. The scenario, 
however, is very different for two and higher dimensions 
as can already be seen from Fig. \%)p where the phase 
space IPR displays a decrease at weak disorder. In this 
case, the momentum space behavior dominates the phase 
space IPR at both, weak and strong disorder. Examples 
of IPRs in real, momentum, and phase space in dimen- 
sions up to d — 3 are given in Fig. 4 of Ref. |2i|. In the 
following, we will distinguish the cases d = 1 and d > 2. 

1. IPR at weak disorder for d = 1 

First, we briefly review the phase space properties of 
the one-dimensional Anderson model, which were already 
discussed in Refs. O and0 in some detail. For W = 0, 
two plane waves at momentum values k and —k are en- 
ergetically degenerate, and there is an ambiguity in the 
choice of the corresponding two basis states. We choose 
symmetric and antisymmetric combinations of the two 
plane waves, in order to obtain real wave functions. The 
solutions in the limit W — > singled out by degener- 
ate perturbation theory contain additional phases which, 
however, do not influence any of the discussed IPRs. 

In the clean case, W = 0, one finds — 1/2 for states 
with non-vanishing momentum. This corresponds to the 



equally weighted contribution of the two momenta k and 
—k. In real space, the sum appearing in the IPR can be 
approximated by an integral and it is suffi cient to con- 
sider as a representative the wave function \/2/L cos(kx) . 
This yields for the IPR in real spaced 




In phase space the Husimi function resolves the two mo- 
menta k and — k which are well separated for energies 
around the band center. While for W — > oo a single 
stripe in phase space leads to an inverse participation 
ratio L -1 / 2 , the two stripes now result in L~ 1 / 2 /2. 

The presence of a disorder potential leads to a cou- 
pling of plane waves with different momenta. In contrast 
to the opposite case of strong disorder, W — > oo, where 
the coupling of the Wannier states occurs only between 
neighboring sites, for weak disorder W — > 0, the cou- 
pling of the plane waves is not restricted to neighboring 
momenta. In fact, the averaged matrix element of the 
disorder potential is independent of the momenta of the 
states involved. Within perturbation theory, however, 
the energy difference of the states comes into play so that 
effectively the coupling to states close in energy (but not 
necessarily in momentum) is dominant. 

Only for the one-dimensional case, the dispersion rela- 
tion E = — 2t cos k implies that real basis states at W = 
which are close in energy E are also close in momen- 
tum k. As a consequence, only states which are close 
in momentum are efficiently coupled by a weak disorder 
potential. Although the perturbative treatment is more 
complicated for weak than for strong disorder, a qualita- 
tive impression of its effect on the phase space properties 
can be obtained in analogy to the case W — > oo by inter- 
changing real and momentum space. Now, because of the 
limited resolution of the Husimi function, the coupling to 
close states in momentum space does not have a signifi- 
cant effect while the large scale modulation in real space 
associated with the coupling affects the Husimi function. 
Therefore, while the value of the IPR in phase space for a 
clean one-dimensional Anderson model at W = is a di- 
rect consequence of the localization in momentum space, 
the corrections for finite W — * are dominated by real 
space effects. 

2. IPR at weak disorder for d > 2 

Generic ballistic states on a d-dimensional cubic lat- 
tice display a 2 d d!-fold energetic degeneracy. The fac- 
tor 2 d arises from the degeneracy between momentum 
vectors with different signs of the components, while 
the factor dl accounts for the number of possible per- 
mutations of a set of d momentum values, provided 
they are all different. For example, in three dimen- 
sions all the 8 different combinations of the signs in 
k = (±fci,±&2, ±&3), and all the 6 permutations of the 
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FIG. 5: Distributions of the IPRs in real space (triangles), 
phase space (circles), and momentum space (squares) in the 
zero disorder limit W — > 0. A generic set of eightfold degen- 
erate states in two dimensions with momentum components 
±57r/24 and ±37r/4 is considered for a system of size L = 48 
where the overlap of the Husimi functions for stripes at fci 
and k2 is negligible. The distributions are obtained by taking 
all eight degenerate states for 10 000 disorder realizations into 
account. 

momentum components {fci, ki, k^} lead to the same en- 
ergy E = — 2i(cos fci + cos k-i + cos fca). Occasionally, the 
degeneracies may even be larger. This is the case when 
the same total energy can be achieved by different sets 
of momentum components. 

While for two degenerate states the limit W — ► leads 
to a universal value for the IPRs, this is no longer true in 
the case of higher degeneracies, where the IPRs depend 
on the disorder realization, even in the limit W — > 0. 
This can already be seen from the existence of different 
types of (real) wave functions. An optimal localization 
in momentum space can be obtained by pairing only two 
plane waves with opposite momenta k and — k, leading 
to 

?/>(x) = (2/i) d / 2 cos(k-x). (14) 

A wide distribution in momentum space is achieved by a 
linear combination of all energetically degenerate states. 
For a generic state this yields 

= ( 0d-2MT d] l/2 E E COS (E^ fc ^V 

(Z d - L > P({k, } ) {„} Vi=i / 

(15) 

where the first sum is to be taken over all permutations 
of the set of (different) momentum components ki while 
the second sum over the set of r\i is to be taken over 
all combinations of factors ±1 with r\\ = +1 kept fixed. 
Only for d — 1 the states ip(x) and <p(x) coincide, hinting 
again at the difference between the weak disorder be- 
havior in one dimension and the subtleties appearing in 
higher dimensions. 

It follows that in two and higher dimensions a non- 
trivial distribution of IPRs already appears in the limit 



W — > 0. Fig. [5] depicts such distributions for a set of 
eightfold degenerate states in two dimensions with mo- 
mentum components ±57r/24 and ±37r/4. The system 
size of L = 48 ensures that the overlap of the Husimi 
functions corresponding to the eight different momentum 
vectors is negligible. 

The two states (1141 and (|15|) help to understand the 
distribution for the momentum space IPR. On the one 
hand, t/>(x) yields the maximum IPR in momentum space 
for real wave functions, Pf. = 1/2. On the other hand, all 
plane waves might be equally weighted as in state </>(x), 
thus leading to an inverse participation ratio of 1/8 in 
momentum space. It turns out that the mixing of the 
plane waves due to a random potential is quite efficient 
thus making the first limit rather improbable. 

For the case presented in Fig. [5] the two momentum 
components fci and &2 are well separated on the scale 
of the phase space resolution in momentum direction 
y/w/L. Consequently, the overlap of the resulting stripes 
in phase space is negligible. Therefore, the distributions 
for the IPR in momentum space and phase space coin- 
cide up to a scaling factor L~ d / 2 which arises from the 
finite width of the Husimi function in the d momentum 
directions. 

In contrast to the behavior in momentum and phase 
space, an equally weighted combination of all energeti- 
cally degenerate plane waves leads to a maximum of the 
IPR in real space. For every pair of different and non- 
zero momentum components k\ and fc.2, which are both 
a multiple of 2-k/L, the IPR for such a state becomes in 
an integral approximation 

Px = ^dx^x) 4 ^^. (16) 

When only one momentum direction contributes, the op- 
posite limit is reached and the IPR in real space becomes 

Px = J dx^(x) 4 - ^ , (17) 

which, up to a factor 1/L, coincides with the result in 
one dimension, cf. Eq. (|13|l . Since the equally weighted 
states described by (fTST) now lead to larger values of the 
IPR, the result in real space is essentially a mirror image 
of the IPRs in momentum and phase space (see Fig. |5J) ■ 

In three dimensions the IPRs for a generic situation of 
48-fold degeneracy may be obtained as well. While L 3 / 2 P 
and Pj. yield values between 1/48 and 1 /2, the IPR in real 
space assumes values between 3/2L 3 and 61/16L 3 . The 
latter values are obtained by a calculation analogous to 
that underlying Eqs. (|16(l and l|17fl . 

We emphasize once more that the states discussed 
above represent the generic states. In addition, there 
exist states where some or all of the momentum compo- 
nents are equal so that the number of degenerate states 
is decreased. On the other hand, in certain cases a given 
total energy can be constructed by different sets of mo- 
mentum components thus giving rise to an increase of 
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the degeneracy. These special states are relevant for a 
detailed description of the complete distribution of IPRs 
for a given system size, which may exhibit a complex 
structure. However, in the limit of large system size, 
the generic states discussed above dominate the distribu- 
tions. 

From the perturbative investigation of the IPR in the 
limits of very strong and very weak disorder, we can con- 
clude that real space properties will only dominate at 
weak disorder strength, if a coupling is induced predom- 
inantly between plane waves close in momentum. As 
long as the distance in momentum is below the momen- 
tum uncertainty y/n/L in phase space, such a coupling 
will become apparent only via the large scale real space 
structure appearing in the Husimi function. Therefore, 
in this scenario which is characteristic for one dimension, 
real space dominates at weak disorder strength and the 
phase space IPR increases with increasing disorder. It is 
only the value of P(W = 0) itself which is determined by 
momentum space properties. 

In two and higher dimensions the picture changes dras- 
tically, because states close in energy are not necessarily 
close in momentum anymore. In this case, the disorder 
potential may scatter plane waves into other momentum 
directions and thus induces a strong mixing in momen- 
tum space. In particular, states of the type (|14() . which 
yield large values for P, are affected by such processes. 
The mixing will thus lead to a decrease of the IPR in 
phase space as a consequence of the dominance of mo- 
mentum space. In contrast, the real space structure will 
appear on relatively small length scales which are typ- 
ically not resolved by the Husimi function. As we will 
see in the following section, this decrease of P in d > 2 
implies the existence of a regime of intermediate disorder 
where phase space is well covered by the Husimi func- 
tion and which can be associated with diffusive behav- 
ior. Furthermore, this scenario opens the possibility of a 
derealization- localization transition. 



VI. SYSTEM SIZE DEPENDENCE 

While the Anderson model in three dimensions exhibits 
a phase transition from dclocalized states at weak poten- 
tial to localized states at strong potential, such a phase 
transition is absent for the Anderson model in two di- 
mensions, where all states are localized in the thermo- 
dynamic limit. Nevertheless, for fixed system size, both 
models show qualitatively the same behavior of the IPR 
in phase space, calling for an analysis of the size depen- 
dence of the IPR in order to check whether or not the 
strong increase of the phase space IPR at the crossover 
between the diffusive and the localized regime evolves 
towards an abrupt jump which would indicate a phase 
transition in the limit L — > oo. 

To this end we plot in Fig. versus the system size 
L, the disorder strengths W at which for the Anderson 
model in one, two and three dimensions maxima and 
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FIG. 6: Values of the disorder strength W for the minima 
(open symbols) and maxima (full symbols) of the inverse par- 
ticipation ratio in phase space, as a function of the system 
size. The circles, triangles and squares are for one, two and 
three dimensions, respectively. 



minima of the phase space IPR occur. The locations 
of minima and maxima are shown as open and full sym- 
bols, respectively. While scaling laws cannot be extracted 
from the data^i one can nevertheless clearly observe the 
direction of the shift of minima and maxima as a func- 
tion of the system size. This reveals an important differ- 
ence between the cases of two and three dimensions. In 
three dimensions, the position of the maximum moves to 
lower disorder values when the system size is increased, 
while the position of the minimum shifts in the oppo- 
site direction. From this trend, one can expect that in 
the limit L — > oo, the positions of the maximum and 
the minimum converge towards the same finite disorder 
value, with the emergence of a non-monotonic step in the 
disorder-dependence of the phase space IPR as a clear 
signature of the Anderson transition in phase space. 

In order to get an estimate of the critical disorder 
strength W c for the Anderson model in d = 3, we de- 
pict in Fig. the change of the phase space IPR as a 
function of the system size for fixed values of the disor- 
der strength. For W < W c , the phase space IPR should 
decrease with increasing L while for W > W c it increases. 
For W = 19.1 (indicated by circles), P will increase with 
the system size. For W = 17 .4 (squares), P decreases for 
the system sizes accessible to us, but one may anticipate 
that the curve will rise for larger system sizes. Such a 
behavior can be observed also for the two-dimensional 
Anderson model. In contrast, for W = 15.8 (triangles), 
one would expect that the curve continues to fall even for 
larger system sizes. This implies a critical value for the 
disorder strength between 15.8 and 17.4. While these 
considerations are not necessarily stringent, the results 
presented in Fig. [7| as well as in Figs. |21 and are per- 
fectly consistent with the known value of W c f=s 16.5 for 
the Anderson transition in the band centeriSi 

In contrast, according to the data depicted in Fig.[S]for 
two dimensions, the positions of the minimum and of the 



9 



1 r 



EN 

m n 1 

0.1 




0.01 



14 16 18 20 22 24 
L 



However, independently of the dimension the increased 
density of states dominates and scattering becomes more 
effective as the system size increases. This corresponds 
to a shrinking of the ballistic regime which can be seen in 
Fig. |2 as a shift of the curves to smaller W with increas- 
ing L. This discussion, however, does not restrict the 
position of the minimum of P as a function of system 
size since the minimum always appears in the diffusive 
regime. Indeed, as Fig. HO shows, with increasing L, the 
position of the minimum clearly shifts towards weaker 
disorder in two dimensions while it shifts to stronger dis- 
order in three dimensions. 



FIG. 7: Dependence of the phase space IPR on the system size 
for the Anderson model in d = 3 at fixed disorder strengths 
W = 6.9, 9.1, 11.0, 12.0, 13.2, 14.5, 15.8 (triangles), 17.4 
(squares), 19.1 (circles), 20.9, 22.9, 25.1, 27.5, 30.2, 36.3, and 
47.9 from the lower to the upper curve. 

maximum IPR both move towards lower disorder values 
when L increases. It seems plausible that they both go 
to zero in the limit of infinite system size, consistent with 
an extension of the localized regime down to infinitesimal 
disorder strength and the absence of a phase transition 
in two dimensions. However, the fact that the overall 
behavior at finite system size in two and three dimensions 
is very similar, hints at the role of d — 2 as a marginal 
dimension in the Anderson model. 

A better insight into the behavior of the phase space 
IPR can be gained by considering the position of the 
maxima, which shift to smaller disorder strength with 
increasing system size, independently of the dimensional- 
ity. For a given system size L, the maximum phase space 
IPR appears in the localized regime at a certain disor- 
der strength W . Now, for a localized state at this fixed 
disorder strength, but at larger system size, the phase 
space IPR becomes independent of the spatial structure 
once the width of the minimal uncertainty state ex- 
ceeds the localization length. In this regime, the phase 
space IPR is dominated by momentum space features. 
Since we know that dPk/dW < we can conclude that 
dP/dW < 0. Therefore, the maximum of P(W) shifts 
to smaller disorder strength when the system size is in- 
creased. Furthermore, from this argument it follows that 
in one and two dimensions the maximum shifts to W = 
in the limit L — > oo if we infer that all states are localized 
in the thermodynamic limit. 

In order to discuss the position of the minimum IPR, 
we now turn to the ballistic regime at weak disorder. The 
coupling between plane waves within first order pertur- 
bation theory depends on two contrary effects. On the 
one hand, the number of plane waves into which scat- 
tering may occur increases with the system size. On 
the other hand, since the disorder potentials at differ- 
ent lattice sites are uncorrelated, the individual coupling 
matrix elements decrease with the size of the system. 



VII. CONCLUSIONS 

In this study of the phase space properties of the An- 
derson model, we have demonstrated the potential imma- 
nent to this approach and, in particular, its advantages 
over approaches based purely on real or momentum space 
properties. In contrast to the latter ones, the phase space 
approach allows to treat real and momentum space on the 
same footing. The well-studied Anderson model has al- 
lowed us to establish an interpretation of the phase space 
IPR which will be useful in cases where no independent 
information is available. 

We found that the crossover between the diffusive and 
the localized regimes is accompanied by an increase of 
the phase space IPR which, in three dimensions, evolves 
to a sharp step in the thermodynamic limit. This is a 
signature of the Anderson metal-insulator transition. 

The jump of the phase space IPR at the Anderson tran- 
sition implies a dramatic reorganization of the Husimi 
distribution from a large spread over phase space to lo- 
calization not only in real space but also in phase space. 
This scenario is not only relevant for the d = 3 Anderson 
model, but corresponds to the very similar one that was 
recently found for the Aubry- Andre model, ^ It is advan- 
tageous to exploit this similarity. The one-dimensional 
Aubry-Andre model allows for a direct visualization of 
the changes in the Husimi function at the metal-insulator 
transition. Furthermore, in numerical treatments of the 
Aubry-Andre model, the system size may be varied by 
more than two orders of magnitude, thus allowing for a 
much more detailed study of the phase transition. 25 

Moreover, by putting together the insights gained from 
phase space into the Aubry-Andre model and the Ander- 
son model at different dimensions, it becomes clear that 
the dimensionality of the model is not the most impor- 
tant parameter for the occurrence of a phase transition. 
Instead, we could identify the disorder-induced coupling 
of plane waves having distant momenta as the relevant 
mechanism for the occurrence of a phase transition. 

It will be interesting to apply these phase space con- 
cepts to interacting systems where the possibility to char- 
acterize individual many-particle states is expected to be 
of great value. Work along these lines is in progress. 
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